Superfluidity of Grain Boundaries in Solid ^He 
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By large-scale quantum Monte Carlo simulations we show that grain boundaries in *He crystals 
are generically superfluid at low temperature, with a transition temperature of the order of ~ 0.5K 
at the melting pressure; non-superfluid grain boundaries are found only for special orientations of 
the grains. We also find that close vicinity to the melting line is not a necessary condition for 
superfluid grain boundaries , and a grain boundary in direct contact with the superfluid liquid at 
the melting curve is found to be mechanically stable and the grain boundary superfluidity observed 
by Sasaki et al. [Science 313, 1098 (2006)] is not just a crack flUed with superfluid. 
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Superfluid grain boundaries (GB) were proposed as 
a plausible scenario [1-3] to explain the effect of non- 
classical rotational inertia (NCRI) in solid ^Hc discovered 
by Kim and Chan [4]. An observation by Rittner and 
Reppy [5] that NCRI signal can be eliminated through 
annealing was the first explicit evidence that crystalline 
defects are of crucial importance. The remarkable direct 
experimental observation of grain boundary superfluidity 
(at the melting point) by Sasaki et al. [6] confirms the 
early theoretical prediction and marks the beginning of a 
new stage in the study of the supersolid phase of helium. 

In this Letter we expand on our previously reported 
preliminary results [3] and show that a grain boundary 
in solid Helium is generically superfluid at low tempera- 
tures. The transition temperature, Tc, is strongly depen- 
dent on the crystallite orientation: while it is typically 
of the order of ^ 0.5 K, grain boundaries with special 
relative orientations of the two grains are found to be in- 
sulating (non-superfluid) . We also obtain strong evidence 
that a grain boundary in contact with the superfluid liq- 
uid, as shown in Fig. 1, is mechanically stable. 

This latter question is important for the interpretation 
of the experimental observation of superflow in crystals 
with grain boundaries [6]. Since these experiments are 
carried out under the conditions of phase coexistence be- 
tween a crystal and a liquid, the observed effect could ei- 
ther be true superflow along a grain boundary, or rather a 
thin liquid-filled crack in a crystal. The answer depends 
on the relationship between the three surface tensions, 
(Ti, (72, and (Tgb (see Fig. 1). Mechanical and energetic 
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FIG. 1: Sketch of the equilibrium configuration of two crystal- 
lites in contact with liquid. The configuration is mechanically 
stable only if the condition (1) is met. Note that surface ten- 
sions (Ti and (72 need not be equal, because of the different 
orientations of the crystalline axes. 

stability of the grain boundary require that 
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If this inequality is not satisfied, the liquid penetrates 
between the two crystallites, and a crack is formed. We 
will show that this configuration is stable. The grain 
boundary superfluidity observed by Sasaki et al. is thus 
not merely the manifestation of a crack in the crystal 
filled with liquid. 

Our Path Integral Monte Carlo (PIMC) simulations 
are based on the continuous-space worm algorithm [7]. 
For spatial imaging, we employ two slightly different 
techniques. The first consists of producing condensate 



FIG. 2: (Color online) Initial setup for the simulation of a 
grain boundary in contact with a liquid. Two truncated pyra- 
mids are placed on top of each other. The basal plane of both 
pyramids is a square with size x Ly = L x L, with L = 24 
[10]. The upper and lower pyramid have different random 
orientations and the height of both pyramids is = L/2. 
The upper facets of the truncated pyramids are squares of 
initial size L/2 x L/2, and form a grain boundary between 
the two crystallites, indicated by the blue square. Liquid fills 
the volume outside the crystallites. Periodic boundaries are 
used in the x and y-directions, while atoms in the 2 = and 
z = L plane, drawn in red, are pinned, to prevent flow along 
the boundary in the z-direction. We refer to the Appendix 
for details such as the exact positions and rotation angles of 
the crystalline lattices within the pyramids. 
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FIG. 3: Condensate map of the two-pyramid sample. The 
map represents (by the density of points) the condensate 
wavefunction in the slice x £ [0.4L, 0.6L] averaged over the 
a;-direction. The initial setup is shown by dashed lines. 



maps, which are maps of the condensate wave function, 
in thin slices of our sample. Within the worm algorithm, 
this is accomplished by recording spatial positions of the 
two open ends of the worms, when they are sufficiently 
far away from each other such that correlations between 
them are negligible [8]. As a result, the density of points 
in the map is proportional to the condensate wave func- 
tion. The second technique is based on the winding-cycle 
maps. Here, we collect statistics of instantaneous particle 
positions, by considering only particles which participate 
in macroscopic exchange cycles characterized by non-zero 
winding numbers, i.e. cycles that directly contribute to 
the superfluid response [9]. 

Stability of the grain-boundary-liquid junction. In or- 
der to check whether a grain boundary is destroyed when 
brought into contact with the liquid, we perform a di- 
rect simulation of two crystallites in contact with liquid 
sketched in Fig. 1. Our simulation setup, shown in Fig. 2, 
consists of two truncated solid pyramids with random 
crystalline orientation placed on top of each other. The 
rest of the volume is filled with liquid. 

Since the goal of this particular simulation is to study 
the stability of a grain boundary we can work at a rel- 
atively high temperature, T=0.8 K, which significantly 
enhances the performance of the Monte Carlo simula- 
tions [11]. The simulation is carried out in the grand 
canonical ensemble, with the chemical potential fixed at 
the phase coexistence point. The equilibrium number of 
atoms in our sample was measured to be about 13660. 



In order to stabilize the solid phase in the system, we 
pin (i.e., do not update) solid atoms in the vicinity of the 
pyramid bases. 

In Fig. 3 we show the condensate map of the sample. 
We sec that the grain boundary between the two crystal- 
lites is a robust quantum object with a thickness of order 
3 (see also Fig. 5) and did not disappear during the sim- 
ulation run. The system has converged to a state where 
the liquid and solid phase coexist. Compared to the ini- 
tial configuration, the shapes of the crystallites (including 
the angles) have noticeably changed. This is not surpris- 
ing since the optimal shape of the crystal-liquid interface 
depends on the particular orientation of the crystallite 
axes with respect to the interface. 

Superfluidity of grain boundaries. Having confirmed 
the stability and superfluidity of grain boundaries in a 
generic sample, we next perform a more systematic study 
of the properties of grain boundaries, extending our pre- 
viously reported simulations [3], where the superfluidity 
had been observed, but only for one particular polycrys- 
talline structure, and grain boundary superfluidity could 
not clearly be distinguished from grain edge superfluidity. 

The simulation of a grain boundary bewteen two trun- 
cated pyramids with basal length L = 24 is computa- 
tionally very expensive since it requires ^ 13660 atoms, 
and is by far the largest simulation of solid Helium per- 
formed to date. Reducing the size of the pyramids is not 
an option, since the typical width of the superfluid grain 
boundary region, estimated from Fig. 3 (see also Fig. 5), 
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FIG. 4: (Color online) Sketch of the initial of two cuboid crys- 
tallites with a total of about 2000 atoms placed on top of each 
other. The basal plane of both cuboids is a square with size 
Lx X Ly = L X L, with L — 12. The upper and lower cuboid 
have different random orientations (see the Appendix for more 
details). The height of both cuboids is Lz/2 — 7, yielding a 
grain boundary in the z = plane (light blue) and a grain 
boundary in the z — Lz/2 plane (dark blue). Periodic bound- 
ary conditions are used along all directions, but motion across 
the xz boundary was suppressed by pinning all atoms at a dis- 
tance smaller than 0.75 from the y = {y = Ly) boundary 
indicated in red. This ensures that any superfluid response 
in the a;-direction is due to the superfluidity of two horizontal 
grain boundaries[12]. An additional grain boundary due to 
periodic boundary conditions arises for each crystallite at the 
X = Q {x = Lx) plane (yellow and orange for the lower and 
the upper cuboid, respectively), but these do not affect the 
superfluid response in the 2:-direction. 



is ~ 3. Hence, in the two-pyramid samples of linear sizes 
significantly smaller than L = 20 30, the grain bound- 
ary is not well defined, and the system is not supposed to 
be metastable. For instance, in a sample with L = 12, we 
found that the grain boundary and the two crystallites 
(apart from the pinned atoms) melt. 

To study the properties of generic grain boundaries 
we use the different sample geometry shown in Fig. 4 
consisting of 2000 atoms . Two cuboids of equal size and 
with random crystalline-lattice orientations are placed on 
top of each other along the z-direction, which creates two 
parallel grain boundaries in the xy-plane due to periodic 
boundary conditions. 

An important result of our simulations is that, as al- 
ready observed in our preliminary studies [3], not all 
grain boundaries in the hep '^He crystal are superfluid. 
Grain boundaries featuring an extra symmetry such as 
stacking faults and special grain boundaries with nicely 
matching angles similar to the one shown in Fig. 6, are 
insulating. 

The generic case, however, are superfluid grain bound- 
aries. Results of a typical simulation at T = 0.25 K and 
n = 0.0287 (melting density) are shown in Fig. 5 . 
The winding-cycle map clearly reveals superfluidity along 
the x-direction. Based on the relatively large system size 




FIG. 5: Phase coherence properties (winding-cycles maps) of 
grain boundaries in the *He sample shown in Fig. 4. Upper 
panel: projection on the xy plane of the data for the upper 
half of the sample containing one of the two grain boundaries. 
Lower panel: (z-shifted) projection of all the data points on 
the xz plane. Note that the x = Q (x = Lx) grain boundary 
(induced by boundary conditions) of one of the two crystal- 
lites turned out to be insulating. This is yet another example 
of insulating grain boundary with special orientation of crys- 
talline axes. 
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FIG. 6: Example of an insulating grain boundary. This grain 
boundary is between two crystallites that differ only by a ro- 
tation about the axis perpendicular to the base plane. Shown 
in the figure is a base plane. The rectangular box is the ac- 
tual simulation cell, the rest being a periodic continuation of 
the sample. Even with this moderate system size, it is clearly 
established that the grain boundary is not superfiuid. De- 
creasing the angle between the two crystallites strengthens 
the insulating character, due to the better match between the 
atoms of different crystallites. In the small-angle limit, this 
grain boundary can be viewed as an array of dislocations; the 
latter are inevitably insulating, given the fact that even at 
moderate angles the grain boundary is not superfiuid. 

utilized in this study, we argue that our observed super- 
fluid signal reflects macroscopic grain boundary superflu- 
idity, and is not an artifact, due for example to the vicin- 
ity of a superfluid ridge — the intersection of the grain 
boundary of interest and the additional grain boundary 
caused by periodic boundary conditions at x = (x = L). 
Incidentally, we note that in the area close to the ridge, 
of size 3, the density of the map is signiflcantly in- 
creased, which is consistent with our general observation 
that nearly all ridges have robust phase-coherence prop- 
erties. 

An interesting observation is that the density in the 
vicinity of superfluid grain boundaries is close to that 
of the crystal, which again confirms that the superflow 
along the grain boundary is not a liquid-filled crack. 

On the basis of a number of simulations similar to the 
one presented in Fig. 5, performed at different pressures, 
temperatures, and crystallite axes orientations, we con- 
clude that generic grain boundaries are superfiuid with 
typical (orientation dependent) transition temperatures 
of about half a Kelvin. The width of the superfiuid grain 
boundary region is ^ 3. We conjecture that the maxi- 
mum possible Tc for grain boundary should be at least 
smaller than the transition temperature of the overpres- 
surized liquid of the same density as the crystal. By sim- 
ulating superfluid properties of the overpressurized liq- 
uid, we found the transition temperature at the density 
n = 0.0287 to be 1.5(1) K. We thus take this value 
as the upper bound for grain-boundary Tc at the melting 
pressure. 

In order to provide a further assessment of the robust- 



ness of our conclusions, we have repeated the same study, 
replacing helium with molecular para- hydrogen, at the 
same low temperatures. In this case, individual particles 
have a mass which is half that of helium atoms, whereas 
the interaction potential is approximately three times 
deeper. No evidence of superfluidity was ever observed in 
this case, in agreement with experimental findings [13]. 

Summarizing, based on a direct quantum Monte Carlo 
simulation of a grain boundary in hep ^He in contact 
with liquid under the conditions of phase coexistence, we 
argue that it is thermodynamically stable against dissolu- 
tion into two crystallites separated by a crack. This lends 
theoretical support to the observation of Sasaki et al. [6] 
of generic superfiuid grain boundaries. We have studied 
superfluid properties of grain boundaries and found that 
special grain boundaries of higher symmetry are insula- 
tors. On the other hand, grain boundaries of a general 
form are found to be superfluid, with typical transition 
temperatures of the order of ^ 0.5 K. 

While these results lend support to an explanation of 
the results of Kim and Chan [4] as grain boundary su- 
perflow, the direct relevance of these findings to the ex- 
periment by Kim and Chan is open, until the presence or 
absence of grain boundaries is investigatred by X-ray or 
neutron scattering experiments on the torsional oscillator 
cells used in these experiments. 
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APPENDIX : PREPARATION OF THE INITIAL 
SAMPLE 

Euler angles 

The rotation from the x, y, z frame to the x'j/'z'-frame 
is completely determined by the three Euler's angles 
(6*, (/), -0). The first rotation is by an angle <^ about 
the z— axis, the second is by an angle 6 G [0,7r] about 
the X— axis, and the third by an angle '0 about the 
z— axis (again). The transformation matrix U, defined 
as r = To + Ur', reads: 

(cV'c0 — c9sips(j) cips(j) + c6s^c(j) s6s^ \ 
—sipc4> — c6cips(j) —stpscj) + cOcipccj) sOctp B A-2) 
S0S(l3 -S0C(j) ce I 
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where c is a shorthand notation for cos and s for sin. 



Construction of the initial sample with the 
pyramids, see Fig. 2 

Based on the crystaUographic definition of the x'y'z'- 
system, the position of the first atom, r' = (0,0,0), un- 
ambiguously defines an infinite hep lattice, of which we 
take only the atoms whose positions are inside the pyra- 
mid. The randomly selected values which specify the 
initial configuration of the lower crystallite are shown in 
Table I. 



The values of the initial parameters of the upper 
cuboid are shown in Table IV. 



parameter 


value 


Xo 


34.946 


yo 


24.507 




0.87081 


e 


1.9978 




4.0671 




2.3174 



parameter 


value 


xo 


34.946 


yo 


24.507 


Zo 


1.0160 


e 


1.9978 


4> 


4.0671 




2.3174 


al parameters for 1 


parameter 


value 


Xo 


6.6684 


yo 


28.603 


Zo 


19.706 


e 


2.1532 




4.7702 




6.0979 



TABLE IV: Initial parameters for the upper cuboid. 



TABLE I: Initial parameters for the lower pyramid. 



The second pyramid is constructed similarly, and then 
mirrored with respect to the z = L/2 plane. The ran- 
domly selected values which specify the initial configura- 
tion of the upper crystallite are shown in Table IL 



parameter 


value 


Xo 


3.9613 


yo 


37.245 




19.288 


e 


0.98923 




1.0833 




1.9885 



TABLE IL Initial parameters for the upper pyramid. 

The initial positions of atoms of the liquid component 
are selected at random. The liquid fills the volume out- 
side the crystallites. 

Construction of the initial sample with the cuboids, 
see Fig. 4 

The initial sample has about 2000 atoms consisting 
of two cuboid crystallites placed on top of each other 
along the z-direction. The crystallites were prepared in 
the same way as the pyramids in Section : we choose 
a random position Tq and random Euler angles (6, (j), ip) 
in order to specify the offset and the orientation of the 
crystallites versus the simulation box. There is no liq- 
uid phase in the intial sample. The values of the initial 
parameters of the lower cuboid are shown in Table III. 
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